The A-Star Puzzle 



Proceedings IAU Symposium No. 224, 2004 © 2004 International Astronomical Union 

J.Zverko, W.W.Weiss, J.Ziznovsky, & S.J. Adelman, eds. DOI: 00.0000/X000000000000000X 



3D-simulation of the Outer Convection-zone 

of an A-star 

Regner Trampedach 1 

1 Research School of Astronomy and Astrophysics, Mt. Stromlo Observatory, Cotter Road, 

Weston ACT 2611, Australia 
email: art@mso.anu.edu.au 



Abstract. The convection code of Nordlund & Stein has been used to evaluate the 3D, radiation- 
coupled convection in a stellar atmosphere with T g = 7300 K, logg = 4.3 and [Fe/H]= 0.0, 
corresponding to a main-sequence A9-star. I will present preliminary comparisons between the 
3D-simulation and a conventional ID stellar structure calculation, and elaborate on the conse- 
quences of the differences. 
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1. Introduction 

From 3 dimensional simulations of convection, it has been known for the last two 
decades that convection grows stronger with increasing effective temperature and de- 
creasing gravity. By stronger, is here meant larger Mach-numbers, larger turbulent- to 
total-pressure ratios and larger convective fluctuations in temperature and density. The 
stronger convection has also been accompanied by increasing departures from ID stellar 
models that fail to predict the extensive overshoot into the high atmosphere, the turbu- 
lent pressure and its effect on the hydrostatic equilibrium, the temperature fluctuations 
and the coupling with the highly non-linear opacity. The latter has the effect of heating 
the layers below the photosphere, thereby expanding the atmosphere, as also done by 
the turbulent pressure. The various ID convection theories/for mulations, e.g ., classical 
mixing- length l|Bohm-Vitenselll958|) . non-lo cal extensions to it llGough1ll977|) or an in- 
dependent formulation based on turbulence l|Canuto 8z Mazzitellilll99^b all have similar 
shortcomings with respect to the simulations. Their predictive power is further limited 
by the free parameters involved. 

Going towards earlier type stars, apart from stronger convection, also means a more 
shallow outer convection zone. This combination is rather unpredictable and is the main 
motivation for the work presented here. Classical predictions call for the outer convection 
zone to disappear close to the transition between A and F stars, but details about where 
and how this transition occurs can only be gained from realistic, 3D simulations, as 
outlined below. 
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2. The simulations 

The simulation presented here was carried out using the cod e of 
|l99ft and is further described in lNordlundl l|l982lb ISteirl l|l989i) and 
(2003). 

The basis of hydrodynamics is the Navier-Stokes equations, and the code employs the 
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conservative or divergence form 
dg 



;H = -V • (qu) (2.1) 



dgu 

dge 



-V • (guu) - VP g + gg (2.2) 
-V • (geu) - P g V • u + e (Q tad + Q visc ) , (2.3) 



where g is the density, P g is the gas pressure, e is the specific internal energy, u is the 
velocity field, g is the gravitational acceleration and Q r ad and <5 V i SC are the radiative and 
viscous heating, respectively, the latter arising from the numerical diffusion applied. 

Equations (|2.1|) - (|2.3|> describe the conservation of mass, momentum and energy, re- 
spectively, with sources and sinks on the right-hand-side. For the convection code the 
equations are preconditioned, by dividing by g, to improve handling of the large density 
contrast between the top and bottom of the simulations. 

The vertical component of the momentum equation, Eq. I|2.2|l . can be written 

F z = - d -^±^ + g 9 , (2.4) 

as we have chosen g to be in the z-direction. With F z = this equation describes 
hydrostatic equilibrium, where the gas pressure and the turbulent pressure, Ptmb = gul, 
provide support against gravity. 

The gas pressure, P g (g,e), and the opacities going into the computation of the radia- 
tive heating, Q ra d, form the atomic physics basis fo r the simulations . The continuous 
opacities are calculat ed from the MARC S-package IjGust afssonl 1 1 9 731) and subsequent 



updates as detailed in lTrampedachl |l997 L the line opacity is in the form of opacity dis- 
tribution functions (ODFs) l|Kuruca ll992). and the equation of state accounts explicitly 
for all ionization st ages of the 15 most abundant elements IjHummer fc Mihalasl ll988: 
iDappen et"allll988|) . 

The present simulation is performed on a 100 x 100 x 82-point grid, has T e s = 7300 K, 
\ogg = 4.3 and [Fe/H]= 0.0, and therefore corresponds to a A9 dwarf on the main- 
sequence. The computational domain is 11.5 Mm on each side, and 13.1 Mm deep, of 
which 1.5 Mm is above the photosphere. 

So far the convection code has only been used for stars that were convective at the 
bottom boundary, so in order to accommodate this simulation, the boundary condition 
was changed and evaluation of radiative heating in optically thick layers was included. 

The simulation was carried out in the plane-parallel appr oximation and includes n o 
rotation or magnetic fields, and has a simple Solar abundance ijAnders fc Grevessell989|) . 
Both thermodynamics and radiative transfer is performed in strict LTE. The velocity- 
field is too large, through-out the simulation domain, to support segregation of elements, 
rendering diffusion and radiative levitation of individual species, comfortably irrelevant. 



3. Comparison with ID models 

In Fig. ^ the temporal and horizontal average of the simulation is compared with a 
corresponding ID stellar model. From the solid black curve in panel a) we see that this 
simulation has an impressive turbulent pressure, making up almost 35% of the total 
pressure about 700km below the photosphere (it is about 13% in Solar simulations). 
The mach-numbers reaching Mach 0.7 as shown with the dashed black line, is equally 
impressive. The convection is much less efficient compared to the Solar case, and invokes 
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both higher velocities, as seen in panel a), and higher super-adiabatic gradient, as seen 
from the black curve in Panel b) of Fig. ^ The two gray curves show the Pt ur b / Ptot-ratio 
and V — V a( j for two ID envelope models using the standard mixing-length formulation 
of convection ijBohm- Vitensdll958|) . The solid line is for a = 1.0 and the dashed line for 
a = 2.0, in an attempt to bracket the behavior of the simulation. One glance at Fig. ^ 
makes it clear that no mixing-length model can reproduce the outer few percent of an A9 
dwarf. The convection zone extends further up into the atmosphere and has much broader 
features than possible with mixing-length models. We also see that overshoot from the 
convection zone supports an appreciable velocity-field, with velocities of 2-3 km s -1 , and 
obvious consequences for spectral line shapes. 

The large Pturb results in smaller temperatures at the same hydrostatic pressure, when 
compared to a ID model. This can be translated into higher pressures and densities on 
an optical depth scale. Consequently the derived population of line-producing states in 
atoms, ions and molecules will be misleading and result in erroneous abundance analysis. 

The mismatch between the simulation and the ID models in the atmosphere is ac- 
tually so profound that no mixing-length model with T c g and g S urf, consistent with the 
simulation, can be found to match the simulation (i.e., P, g and T) at the bottom. This 
means that stellar structure and evolution calculations are misplaced in the HR-diagram, 
with implications for age-determinations of globular clusters and our general knowledge 
of the interior and evolution of A-stars. 

4. Temperature inversion 

This simulation, which is the hottest performed with this code, displays some large 
and persistent temperature inversions, 1.5-2.5 Mm below the photosphere. They have 
an amplitude of up to 8 000K and are about 0.5 Mm deep. During the meeting, Noels 
suggested that there might be a correlation between the turbulent pressure and the low 
temperatures; keeping P to t and g constant to ensure hydrostatic equilibrium, a large 
Pturb would force a low temperature. This is, however not observed in the simulations. 
There is no correlation between the temperature inversion and Pturb ~~ nor are there cor- 
relations with vertical or horizontal velocities, or the density contrast. There is, however, 
a strong correlation with the vertical force (see Ea. l2.4fl . as the force on the plasma in the 
temperature inversion is always outwards (the — z-direction). Furthermore, the vorticity 
is clearly skewed towards larger values in the inversion layer, especially when compared 
to the upflow, but also with respect to the down-drafts. 

Fig. El shows a vertical snapshot in temperature, displaying a typical temperature in- 
version. The arrows indicate the velocity field, with the maximum length corresponding 
to 2.9 km s -1 . The temperature inversion extends some 6 Mm across the left side of the 
plot, at a height of 2 Mm. It is connected to the downdraft at the edge, but it never 
connects with the downdraft to the right of the center. The velocity field there is con- 
nected, though, but the colder plasma from the inversion is compressed on the way to 
the downdraft, and therefore heated. The velocity field is diverging in the z-direction, 
neutral in the x-direction and convergent in the y-direction, perpendicular to the plane 
of the figure - the net-flow is convergent. The features observed in this simulation are 
cylindrical in the horizontal direction, with width and depth being similar and the length 
being 5-10 times larger. 

The temperature inversions seem to develop when the local photosphere has subsided 
to about 2 Mm below the average photosphere, at the edge of a granule. The layers above, 
closer to the height of the average photosphere, then heats up, leaving a cooler area in 
between - the temperature inversion. In white-light surface intensity, this now looks like 
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Figure 1. Comparisons between the simulation and two ID, mixing-length models. In both 
panels the black solid line shows the temporal- and horizontally averaged simulation, and the 
gray solid line shows the same quantity for a ID model with a — 1.0 and the dashed gray curve 
is for a = 2.0. Panel a) shows the turbulent- to total-pressure ratios, and the dashed black line 
shows the RMS Mach-numbers. Panel b) shows the super-adiabatic gradient, with the dotted 
zero-line aiding the location of the borders of the convection zones. 



the edge of a normal granule. The inversion immediately starts eroding from the newly 
heated region on top moving down, and presumably from heating by the surroundings. 
The temperature profile after this sequence, looks like that of a downdraft, and from 
the surface, what looked like the edge of a normal granule, collapse as the temperature 
inversion below disappear. The whole cycle takes about 5 min. 

The reason for this behavior is still under investigation, but it might be connected 
to the local minimum in V s , as seen in Fig. ^ around log 10 P g ~ 4.9 corresponding to 
z ~ 1.9 Mm. 



5. Summary 

A realistic 3D simulation of convection in the surface layers of a A9 dwarf reveals pro- 
found differences between that and a conventional ID stellar structure model, as indicated 
in Sect. |21 In general, the simulations have much smoother and broader convection-zone 
features, compared to mixing-length models. 

We are also treated to a new phenomena, as the simulations display repeated local 
temperature inversions just below the photosphere (c/., Sect.^J. The mechanism respon- 
sible for these large inversions, have not been uncovered yet, but it being an effect of 
large fluctuations in the turbulent pressure, has been ruled out. 
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Figure 2. A vertical slice of a snapshot showing a temperature inversion. The color-scale 
shows temperature, and the arrows indicate the velocity-field. 

This is still work in progress, and in the future, this and other simulations will be used 
for evaluation of spectral lines, limb-darkening, broad-band colors, granulation-spectra 
and p-mode spectra. 
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Discussion 

COWLEY: Did you get a higher contrast between hot and cold regions than in simulations 
of cooler stars? At one time I though the ionization imbalance that we observe in some 
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cool Ap stars might be explained by such temperature inhomogeneities. But I did some 
test calculations and found this could not explain the ionization - at least not with the 
numerical models I had. 

Trampedach: The RMS temperature fluctuations at r = 1 on a Rosseland optical depth 
scale, are about twice as large for the A9 simulation (~ 610 K) as for the Solar simulation 
(~ 300 K). For the latter, the RMS fluctuations stay below 750 K, peaking at r = 10, 
whereas in the A9 simulation they grow almost exponentially to reach 5000 K at r = 100. 

I don't know what kind of numerical models you have used, but a decisive factor in the 
emergent spectra from the convection simulations is the asymmetry of the distributions of 
temperatures, densities and velocities, as well as their correlations. These characteristics 
are hard to predict without realistic simulations. 

PiSKUNOV: What is the height of the temperature inversion and the filling factor of the 
effect? Do you see significant changes of the filling factor? 

TRAMPEDACH: The temperature inversions are about 0.5-1.0 Mm deep and peak at a 
depths of 1.5-2.5 Mm. They involves from 2% to 8% of the horizontal area, with an 
average of about 5%. 

Noels: Do you find a high turbulent pressure in the region of temperature inversion? If 
so, I would suggest to investigate this point in relation to the origin of the temperature 
inversion. 

TRAMPEDACH: There is no obvious correlation between the temperature inversion and 
the local turbulent pressure. I would refer to Sect.0]for the present extent of my analysis. 

Grevesse: Micro-turbulence and macro-turbulence needed with ID stellar photosphere 
models are not needed anymore with 3D models. Why is it not the same when 3D models 
are used, instead of ID models, for A stars? 

Trampedach: I don't think that has changed, going from solar-like stars, to A stars. 
I think what you are alluding to, is the lack of reversal of bisector shape in the A 
star simulations presented by Freytag (elsewhere in these proceedings). I don't know 
the reason for this and the only major thing missing from those simulations is line- 
blanketing. Whether that will make the difference is unclear and will have to await 
further simulations. It is important to understand that theoretical bisectors from ID 
models are straight, no matter how much micro- or macro-turbulence is applied. The 
shape of bisectors is a higher order problem, that cannot be addressed with ID models. 



